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Abstract 

We study numerically the evolution of wavepackets in quasi one-dimensional 
random systems described by a tight-binding Hamiltonian with long-range 
random interactions. Results are presented for the scaling properties of the 
width of packets in three time regimes: ballistic, diffusive and localized. Par- 
ticular attention is given to the fluctuations of packet widths in both the 
diffusive and localized regime. Scaling properties of the steady-state distribu- 
tion are also analyzed and compared with a theoretical expression borrowed 
from the one-dimensional Anderson theory. Analogies and differences with 
the kicked rotator model and the one-dimensional localization are discussed. 
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I. INTRODUCTION 



The main approach to a statistical description of the spectra in complex quantum systems 
originates from the pioneering work of Wigner ffl], who conjectured that random matrices 
could represent the simplest meaningful model for studying heavy nuclei. Currently, ran- 
dom matrix theory (RMT) has become a very effective tool in a large variety of physical 
applications. Until recently the matrices in this theory were assumed to be homogeneous, 
i.e. all matrix elements were taken to have identical statistics. This simplification is mainly 
dictated by mathematical reasons since the corresponding ensembles of random matrices are 
rotational invariant, a property that simplifies the theoretical analysis. 

For a long time, the RMT had no concrete physical basis, in the sense that conditions 
for its applicability were not specified. It was believed that the systems under consideration 
had to be extremely complex in order to have a good agreement with the predictions of the 
RMT [0-0]. The situation has changed with the progress of the so-called quantum chaos 
theory which deals with dynamical Hamiltonian systems exhibiting chaotic motion in the 
classical limit. One of the main results of this theory is that in the extreme case, when 
classical motion is strongly chaotic and no influence of quantum localization is taken into 
account, statistical properties of both spectra and eigenfunctions are well described by the 
RMT. This statement has been thoroughly studied and confirmed both for autonomous 
systems like chaotic billiards [|J and for time- dependent models like the Kicked Rotator 
model (KR) |||6j and the Kicked Tops 0. Moreover, there are many physical examples 
where "complexity" of a quantum system is not maximal, but nevertheless a statistical 
description applies pretty well. 

To describe the consequences of quantum localization in the presence of strong classical 
chaos, a new type of random matrices has been introduced when studying the KR [§]. The 
distinctive feature of these matrices is their band-like structure, which is related to the finite 
range of interactions in a given basis. In a sense, the ensembles of Band Random Matrices 
(BRM) can be regarded as an extension of the conventional random ensembles, since the 
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latter are recovered by just setting the band size b equal to the matrix size N. Currently, the 
interest for BRM raised significantly due to their close relationship with quasi-lD models 
with random potentials ||, or ID systems with long-range hoppings between neighboring 
sites. In the ID interpretation, the band size b corresponds to the hopping range while in 
the quasi-lD case, it is associated with the number of transverse channels for electron wave 
propagation along a thin wire. Recent numerical and analytical studies of BRM (see PJlO 



and references therein), led to numerous results regarding the structure of eigenstates. In 
particular, the localization length has been shown to depend only on the scaling parameter 
b 2 /N, so as the statistical properties of the eigenstates that are directly related to the 
fluctuation properties of the conductance. 

However, much less is known about the evolution of wave packets in models described 
by BRM. One should note that even in the "simple" case of Anderson-type models in ID 
only the short and the long time scales, corresponding to ballistic spread and saturation of 
the packet width, respectively, have been successfully studied [PTf| . In quasi-lD (or ID with 
long-range hoppings) models, the picture is both more complicated and more interesting 
with respect to the Anderson case. Indeed, while the diffusion time scale is absent in ID 
models of Anderson type, since the mean free path is of the order of the localization length, 
classical diffusion alters the ballistic spread in quasi ID systems, before being eventually 
suppressed by localization effects. 

In dynamical systems, the influence of strong localization on classical diffusion in mo- 
mentum (or energy) space has been studied in detail in the framework of the KR. It was 
found that the time scale of the wave-packet spreading that is analogous to classical diffu- 
sion is much longer than the logarithmic time scale over which the complete correspondence 
between classical and quantum description holds 0,^,|1^,|T^] . The entire diffusive process, 



including the final saturation proved to exhibit remarkable scaling properties ||[T3| . In par- 
ticular, the diffusive time scale is proportional to the localization length of those eigenstates 
which are involved in the dynamics. In view of the analogies existing between dynamical 
and Anderson-type localization, the results obtained in the study of wave packet evolution 



in the KR will represent the touchstone for the present investigation of packet dynamics 
in quasi-lD disordered models. However, because of the existence of basic differences as 



ID geometry, with particular attention to their width and fluctuations on different time 
scales. We hope that our detailed numerical investigations, providing direct information 
on the dynamical properties of strong localization, will furnish also some insight for future 
analytical studies. 

The paper is structured in the following way. In Sec. 2, the model is introduced and 
discussed together with the main properties of spectra and eigenstates of BRM. In Sec. 3 
both ballistic and diffusive time scales are analyzed and the scaling properties of packets in 
terms of the band-width b are established. The effect of noise in the destruction of coherence 
is also briefly discussed. Section 3 terminates with the results for the suppression of classical 
diffusion due to the localization of eigenstates. In section 4, we focus on the problem of 
fluctuations of the shape of packets both for the diffusive and relaxation time scales. In 
section 5, a phenomenological description of the asymptotic shape of the packets is given, 
based on results for the ID Anderson model. Moreover, fluctuations along the asymptotic 
profile of packets are studied. The conclusions are summarized in Sec. 6. 



Since its birth, random matrix theory has been mainly dealing with statistical properties 
of "full" random matrices, for which all off-diagonal matrix elements are independent and 
distributed according to the same law. In physical applications this implies that interactions 
are assumed to be so strong and complex that no other parameter, apart from the symmetry 




presently clear to what extent the similarity between 



dynamical and disordered systems can be pushed forward. 



In this paper we extend a previous study jljj of initially <5-like wave packets in a quasi- 



II. BAND RANDOM MATRICES: MAIN PROPERTIES 



A. Definition and applications 



4 



of matrices, is to be taken into account. As a result, such matrices are associated with the 
extreme case of maximal chaos which is known to appear in various physical systems such 
as heavy nuclei, atoms, metallic clusters, etc. Furthermore, full random matrices represent 
a good model for the description of local statistical properties of spectra and eigenstates in 
some range of the energy spectrum, typically, in the semiclassical region. 

On the other hand, the conventional RMT is both unable to describe important phenom- 
ena such as localization of eigenstates, and to characterize the spectra of physical systems 
influenced by strong localization effects. For this reason, much attention has been recently 
paid to the so-called Band Random Matrices which are characterized by the free param- 
eter b defining the effective band width of a Hamiltonian. Such random matrices with 
elements decaying away from the main diagonal, appear to provide more realistic models 
for the Hamiltonian of "complex" quantum systems (see, e.g., ||,|10| and references therein). 
The simplest type of BRM is given by matrices H nm with zero elements outside the band 
(|n — m\ > b), while inside the band (\n — m\ < b), matrix elements are assumed to be 
independent and distributed according to a Gaussian law, 

V(H nm ) = X -= ex V (-HlJ2a 2 nm ), a 2 nm =< H 2 nrn >= ^(1 + S nm ) (1) 
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2^]- — nmi — nmn ~ nm — nm ^ 



where 5 nm is the Kronecker symbol, and Oq = 2 implying that the variance of the off-diagonal 
matrix elements is equal to 1. The BRM ensemble can be regarded as a generalization of 
the standard Gaussian orthogonal ensemble (GOE) as the former reduces to the latter 
when b = N. Analogous generalizations can be introduced for the unitary and symplectic 
ensembles of band random matrices . 

The limit case b = corresponds to diagonal matrices, while 6=1 corresponds to 
tridiagonal Hamiltonians with both diagonal and off-diagonal disorder. The latter case is 
well known in the physics of disordered media; the main properties of such matrices are 
relatively well understood. The general case of BRM where the variance of the matrix 



elements decreases with the distance from the main diagonal, introduced in NM, is also 



amenable to an analytic treatment. They are not, strictly speaking band matrices, but an 



effective band size b can be denned from the shape of the envelope |19 . 



In what follows, we consider large values of b 3> 1 and assume that iV > b 2 . The first 
condition implies a large number of nearby states coupled by the interaction. The second 
condition allows to neglect finite size corrections arising from the finiteness of the samples. 

Considerable interest towards the ensemble of BRM was stimulated by the investigation 
of the quantum behavior of periodically driven Hamiltonian systems. A paradigmatic system 
in this class is the so-called Kicked Rotator model. Indeed, the unitary matrix U yielding the 
time evolution between two consecutive kicks has, in the angular momentum representation, 
a band-like structure with an effective band-size approximately equal to the strength k of the 
kicks. Outside the band, the matrix elements of U decrease extremely fast, while inside the 
band they can be treated as pseudo-random entries if the corresponding classical evolution 
is chaotic 0. As a consequence, both spectra and eigenstates of BRM of the type ([[]) 
are expected to have statistical properties similar to those of the KR. A number of data 
substantiate this belief . 



BRM play also an important role in the understanding of quasi-lD disordered media. 
The connection has been established through the nonlinear sigma model which, on the one 
hand, is known to provide an excellent description of quasi-lD systems and, on the other, 
has been rigorously shown to be reducible to the BRM model Such relationships gave 



a boost to the investigations of statistical properties of eigenstates and eigenvalues of BRM 
as they also allow a better understanding of the properties of the KR, as well as of quasi- ID 
and ID models with long-range random interactions. 

B. Density of States and Structure of Eigenstates 



As was first numerically shown in [z(J and later analytically proved in ||18| , the density 
of states p(E) for infinite BRM (N — > oo) obeys the semicircle law 

1 



p(E) = — -VSbv 2 - E 2 ; \E\ <R = vV8b (2) 
Aixbv 1 
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with p(E) = for \E\ > Rq . The parameter v is just the standard deviation of the 
distribution of off-diagonal elements, v 2 = (H 2 m ) = (Tq/2; it does not influence the statistical 
properties of the spectra as well as the structure of the eigenstates, since it can be scaled 
out. For b = N, expression (Q) reduces to the well known Wigner semicircle law derived in 
conventional RMT. 

In infinite BRM, all eigenstates <fE(n) are known to be exponentially localized [|9jl0,20 



\<p E (n)\ ~ exp I- j j ; n^±oc (3) 

where no is the "center" of an eigenfunction in the basis in which the random matrices 
have been defined. The quantity loa(E) is the localization length defined as the inverse 
of the asymptotic spatial decay-rate of the amplitude of the corresponding eigenfunction. 
Numerically, loo{E) can be determined by implementing the transfer matrix method. It is 
important to recall that the localization length describes the decay of the eigenfunction 
only in the tail and not in the central region of size ~ b 2 . This region is characterized by 
an effective number I ~ of "principal components" which are usually defined in terms of 
the inverse participation ratio and of the so-called entropic localization length Ih (see 
for details). 

In finite samples, one more parameter comes into play, the rank N of the matrices. 
In such a case, all relevant properties of spectra and eigenstates are parametrized by the 
ratio A = b 2 /N PJlQl. Upon changing A, one can accurately follow the transition from 



the completely localized (A <C 1) to the delocalized (A 3> 1) regime. Numerous studies of 
BRM allowed to unravel the dependence of the statistical properties of eigenfunctions on 
this scaling parameter. In particular, in finite bases of size N, all eigenstates are extended if 
A>1, and all their properties are very similar to those for the standard RMT. Our interest 
here, is devoted to the opposite limit A 1 of very localized eigenstates where finite-size 
effects can be neglected. Based on the results obtained for the KR PJ2"0[|, it was predicted 
that the localization length is proportional to b 2 . A rigorous analysis [TE| has confirmed 



this prediction and established the dependence ~ p b of the localization length on the 



energy. In the case of BRM with a general envelope function a(k) for matrix elements H nm , 
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a(k) k 2 

UE) = 2[1 - {^) 2 ]B ■ B= k -^ ; k = n-m (4) 

Ro [ £ a(k)} 2 

k=—oo 

where B is the second moment of the function a(k). In case ([TJ), i.e. for a sharp band of size 
b, one obtains B = b 2 /3 and thus 

U£) = jj[l-(^) 2 ]& 2 • (5) 

At variance with the KR, the localization length of the eigenstates of BRM depends on the 
energy E. 

C. Numerical procedure 

Although considerable progress has been made in the description of the eigenstate struc- 
ture (see |5,ni] and references therein), the evolution of wave packets is still poorly under- 



stood even in the limit of infinite samples. The mathematical model we consider below is 
the time-dependent Schrodinger equation on a ID lattice, 

dr (t) n+b 

^jr= E H nmCm (t), (6) 

' ' m=n—b 

where c n (t) is the probability amplitude for an electron to be at site n and H nm is a symmetric 
BRM of the type (|]). Equation (^) has been integrated numerically using a finite-time-step 
(dt ~ 10~ 3 — 10~ 4 ) fourth order Runge-Kutta algorithm on a self-expanding lattice in order to 



eliminate finite-size effects ||21|| . Whenever the probability of finding the particle at the edges 
of the chain exceeded 10 -15 , 106 new sites were added to each edge. The initial condition 
was taken to be a 5-like state located in the middle of the chain, i.e. c n (t = 0) = 5 n fi. At 
each time step, the normalization condition for the total probability, J2 n l c n(^)| 2 = 1 ; was 
checked observing fluctuations smaller than 10~ 4 . 

A further check of the accuracy of our calculations has been performed by reversing the 
time-axis direction after 2000 time units (for b = 10). The difference between the initial 
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probability distribution and that obtained after integrating for 2000+2000 units was found 
to be less than 10~ 13 . In all the cases discussed below, a large number of disorder realizations 
has been considered (more than 150) in order to get rid of sample-to-sample fluctuations. 



III. DIFFUSION OF WAVE PACKETS 

The time-evolution of a quantum wave packet in the lattice is naturally described by the 
mean square displacement, 

M(M) = («(*)> = ^> 2 M*)I 2 ) , (7) 

where (...) stands for the average over different realizations of the frozen disorder H nm . The 
time dependence of M(b, t) provides a qualitative description of the dynamical regime: a 
power-law evolution, M(t) ~ t u , with v < 1 corresponds to a sub-diffusive behaviour (hinting 
at a possible, eventual localization), v — 1 corresponds to ordinary diffusion, while v > 1 to 
super- diffusion (y — 2 characterizes ballistic motion). 

As mentioned above, BRM can be regarded as a good model for dynamical quantum 
systems such as the KR in the region of strong classical chaos. In the classical limit, this 
model exhibits an unbounded diffusion in angular momentum space if the kick strength 
exceeds some critical value. It was discovered that even in the deep semi-classical domain, 
quantum effects can suppress classical diffusion || giving rise to a phenomenon that is 
closely related to Anderson localization of a quantum particle in random potentials |22]-|24 



This effect of "dynamical localization" was claimed to be experimentally observed in the 
ionization of hydrogen subject to a monochromatic field |25j. A formal connection with a 



ID tight-binding model has been found in [pfj , thus reviving a general interest for localization 
in one-dimensional systems. 

We have investigated the behaviour of M(b, t) by numerically integrating Eq. (||) for 
different values of b. The results of our analysis have been compared with both theoretical 
predictions of the ID Anderson model and numerical findings for the KR. 
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A. Ballistic time scale 



The essential difference between periodic and disordered quantum lattice structures 
mostly lies in the localization properties of their electronic states. In the periodic case, 
all the states are perfectly extended Bloch waves, while in strongly disordered samples, the 
states are asymptotically localized in time because of quantum interference effects. How- 
ever, even in the latter case, there exists a ballistic regime, occurring on time scales of the 
order of the elastic scattering time i.e the time for an electron to move by an amount 
corresponding to the mean free path l m . In quasi-lD systems, l m is known to be equal to 
the number of transverse channels. 

In order to investigate the scaling behavior of the packet size with b, we have numerically 
integrated Eq. @ for very short times ~ tf, and different values of b in the range b = 20-=-45. 
The behaviour of the mean-square displacement M(b, t) is reported in Fig. 1 with the scaling 
assumption 

M(b,t) =b 2 M(tVb) . (8) 

The very good data collapse confirms the scaling Anzatz. Eq. (§) can be understood by 
estimating the ballistic time scale t&. Let us start by noticing that the leading contribution 
to the wavepacket spreading over short time scales comes from the 2b sites which are directly 
coupled with the site n = where the wave packet is concentrated at time t = 0. By 
evaluating the Schrodinger equation and thereby determining M(b,t), one obtains 

M(M)= £ n 2 \H nfi \H 2 « bH 2 . (9) 

n=— b 

The above type of evolution terminates when the average width a/M of the packet becomes 
of the same order as the band-size b, so that farther sites come into play. By substituting 
back in Eq. (||), one finds that the ballistic spread occurs on the time scale 

t<t b n—. (10) 
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Accordingly, the ballistic time scale shrinks to zero for increasing the interaction range b. 
Notice that the ballistic regime is entirely new with respect to the analogous problem in the 
KR where, at small times, an exponentially fast spread of the packet takes place. 



B. Diffusive time scale 

In ID Anderson-type models, wave packet saturation starts immediately beyond the 
ballistic time scale, since the mean free path l m and the localization length are of the 
same order (l^ « 2l m ). Accordingly, no intermediate diffusive regime can be observed. In 
the case of large band-size, b ^> 1, the localization length ~ b 2 is much larger than the 
mean free path l m ~ b, so that a diffusive time scale appears. In order to estimate tu, 
we shall follow the scaling arguments developed from the theory of quantum chaos where 
they have been successfully introduced to explain the well known "quantum suppression of 



classical diffusion" in the KR [|T^]. 



The first crucial observation is that the eigenfunctions fE{n) are exponentially localized 
in the standard basis for all energy values E in the spectrum (consequently, there is a pure- 
point spectrum p{E)). Let us proceed by noticing that an initial state c n (t = 0) excites an 
effective, finite number N e ff of eigenstates with corresponding energies Ei. Accordingly, the 
spectrum of those eigenstates participating to the evolution of the packet is characterized 
by a mean level spacing A e // ~ Ro/N e ff where = 8bv 2 is the radius of the semicircle, 
see Eq.(|). 

Therefore, for times t < 1/A e ff, the packet evolution does not "feel" the discreteness of 
the spectrum, which is resolved over longer time scales and eventually leads to localization. A 
typical evolution is reported in Fig. 2, where the dependence of the mean square displacement 
M(b, t) for b = 12 is shown on a very large time scale. However, as long as tf, < t < tp, with 

t D ^N eff /R (11) 

the motion of the particle is analogous to the standard (classical) diffusion. This means that 
M fa Dt, where D is the diffusion constant and M can be interpreted as the square of the 
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number of effectively excited, unperturbed states. This regime is clearly seen in the inset of 
Fig. 2, where the evolution of M is reported in a doubly logarithmic plot. 

The quantity M reaches its maximal value M max at t ~ to- The value M^ 2 X is °f t ne 
same order as the total number N e jf of eigenstates that participate to the packet evolution. 
Let us finally notice that the packet-width is asymptotically of the order of the localization 
length of the eigenf unctions, i.e. N e ff ~ ~ b 2 (the energy dependence is here irrelevant 
and can be dropped). Accordingly, the following scaling relations hold 

t D ~loo/Ro D-l^Ro , (12) 

where Rq is the width of the spectrum. The second estimate in ( |T2"D corresponds to the well 
known relation between the localization length and the diffusion coefficient in the theory of 
disordered solids, ~ npD where p is the density of states (see, for example, 0). 

Eq. (|L2"|) seems to suggest that the proper scaling of the time axis is 6 3//2 . However, as 
we shall see in the next section, there is convincing evidence that best scaling Ansatz, in 
comparison with the direct integration of the Schrodinger equation (^), is t/b 2 . Whether the 
discrepancy is to be attributed to some weakness of the above arguments or to an improper 
choice of the time units in the numerical procedure it is not clear (in fact, everything works 
perfectly as if we had to refer to dimensionless time units in Eq. @ normalized to <jq - see 
Eq. ([[]) - rather than to the width of the spectrum R ). 

C. Diffusion suppression and scaling properties 

One of the most important peculiarities of the time evolution of wave packets in ID 
and quasi-lD random potentials is the saturation of the width M for t — > oo. In analogy 
with the evolution of wave packets in the KR 0] and from the localization properties of the 
stationary problem (the localization length grows as b 2 ), one expects that for b ^> 1 the 
limiting value M^b) grows as 6 4 . In order to confirm this prediction, we have performed 
detailed numerical experiments in the range b = 4 4- 12. The asymptotic value M 0O (fo) has 
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been accurately determined by averaging M(b,t) over a long time after an initial transient. 
From our data we have found that the dependence of on b is slightly slower than 
expected, ~ b a with a ~ 3.87 ± 0.02. This anomalous behaviour is presumably to 
be attributed to the presence of finite band-size corrections which are not negligible in the 
range of b values that has been numerically investigated (b < 12). Our results are reported 
in Fig. 3, where M(b,t) and t are divided by the asymptotic value and its square root, 
respectively. While the scaling Ansatz for M follows straightforwardly from the detailed 
knowledge of the localization properties, the rescaling of time axis is mainly justified a 
posteriori from the resulting good data collapse that can be observed in Fig. 3. 

Anyway, the nontrivial evolution during the late stages of the diffusive regime confirmed 
by a direct investigation. Since generic properties of eigenstates in quasi-lD models have 
been found to be similar to those of strictly ID disordered models 0, it is natural to expect 
that the similarity extends to the dynamics of wave packets as well. However, even in 1D- 
geometry, the analytical treatment is very difficult. Analytical results are available only in 
the two opposite regimes, t <C 1 and t ^> t D . For example, the asymptotic dependence of 
the mean squared displacement (u(t)) = (x 2 (t)) of packets in the long-time limit is given by 
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(u(t)) ~ a l 2 m ( 1 - ln |/2tl m) ) ; t>>tm (13) 

where l m and t m are the mean free path and the corresponding time between consecutive 
back scattering processes. This estimate is based on the expression for the quantum diffu- 
sion coefficient obtained in p4 |. The logarithmic singularity in (pT3|) follows from resonant 



transitions occurring between pairs of the so-called Mott states |27|]28|| . Such states have a 
peculiar structure characterized by two humps lying at a distance much larger than their 
effective width. Since Mott states appear in pairs, the corresponding energies are very close 
to each other and this results in a resonant tunneling over large distances. The influence of 
Mott states on electronic properties of disordered ID models has been studied in , where 
the clustering of energy levels was discovered and attributed to these states. 

13 



The effect of these states has been included in the study of the long-time behavior of 



wave packets in the KR , where an expression similar to Eq. flT3| ) has been introduced to 
describe the evolution of the mean square displacement M(t) in the momentum representa- 
tion, 

& m^Mt/m). t>>2D (14) 

dt (t/(2D)f V ; 

where D is the classical diffusion coefficient. Numerical data seems to confirm the expec- 
tations (Eq. flliD ), although the presence of very large fluctuations prevent to draw a con- 
vincing conclusion. Despite the better statistics of our data, the presence of the logarithmic 
correction cannot be definitely assessed in BRM too. 

Another approach to the problem of quantum diffusion in the presence of strong localiza- 
tion has been recently suggested in (see further developments in fl4]|): it is essentially 
based on a phenomenological diffusion equation for the Green function, which takes into 
account backward scattering. At large times, the relaxation rate is given by |T3| 



(t/2D) 2 K J 

which differs from Eq. ( |T4"D by a further logarithmic factor. At the moment, all available 
numerical data for the KR do not allow to draw a final conclusion in favor of either expression. 
In any case, let us again remark that quantum localization in the KR is of dynamical nature 
(there is no randomness in the model), so that it is not clear to what extent it is similar to 
the localization of Anderson type. 

Instead of focussing on the question of the exact asymptotic dependence {t — > oo) of the 
mean square displacement M(b,t), it is, for the time being, more useful to limit ourselves 
to provide an effective description of the wide time region that includes also the crossover 
from classical diffusion to complete saturation. In the absence of any theory, we make use 
of the phenomenological expression suggested in |32| (see also B) 



MtM^wfi-^d— ), (16) 
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where M m , to and (3 are the three independent parameters to be determined. The first one 
is obviously obtained from the asymptotic evolution, while the short-time classical diffusion 
(here, we neglect the ballistic time scale which is indeed negligible for large b) provides a 
further constraint to be fulfilled. In fact, for t -C to, Eq. (|IIJ) reduces to 

M (b,t) = = Dt , (17) 

tD 

which allows expressing to in terms of the last unknown j3 which can be determined by 
fitting the global behaviour of M(b,t). 

The main idea behind the phenomenological expression (|i6|) is the repulsion of the energy 
levels participating to the evolution of the packet. As it was previously discussed, the 
diffusion rate is proportional to the mean spectral density and it remains unchanged for 
t < to, according to the uncertainty principle. However, for t > tp it decreases, since 
the only eigenstates that continue to contribute ("operative eigenstates" ) are those whose 
energy level-spacing s satisfies the relation s < t^/t. The relative number of such spacings 
(hence, the relative diffusion rate) is given by the spacing distribution p(s) for the operative 
eigenstates, 

where t ^> tjj is assumed. The above time dependence for large times is the core of the 
phenomenological expression (|T6|) (see details in ||). 

According to the above relation, the parameter (3 characterizes an effective repulsion 
between those eigenstates which are excited by the initial wave packet. It is clear that 
these eigenstates strongly overlap. As a result, the value of (3 can be expected to be quite 
close to 1. Although these arguments are no longer valid for very long times, when the 
level clustering fll6| , |33| due to the influence of Mott states becomes important, Eq. ([H]) can 
still provide a sufficiently good description of the packet dynamics. Numerical experiments 
done for the KR |6|,|32| yield quite a small value of (3 {(3 ~ 0.3). This result, which is 
somehow contradictory with other studies (see the discussion in ) , is probably due to the 
insufficiently long times considered in the simulations. 
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Our detailed numerical experiments with BRM, performed on a much longer time scale 
and with high statistics, reveal quite a good correspondence with the scaling dependence 
(|T6|) (see Fig. 3). The best fit of Eq. (|TJ) gives the following values: ft pa 0.8 and t D ~ 2b 2 . 
Since the values of the band size b are not very large, the limiting value has been 
purposely rescaled to the same level for the different 6-values. By neglecting the residual 
weak deviations from a perfect scaling, we obtain pa 1.9b 2 . The most important point 
of the above analysis is that the value of the repulsion parameter (5 is quite close to 1. 
This means that even for very large times close to the relaxation, the approximate power 
dependence l/£ ' 8 for the difference AM = M (X) (6) — M(b,t) mimics the correct dependence 
AM ~ ln(t)/£ (see Eq. fll"3|)). As a result, one can treat the scaling dependence ( |16D as a 
good description of both classical diffusion and its suppression due to strong localization of 
eigenstates. 

The asymptotic localization of the wave packet is entirely a consequence of the frozen 
character of the disorder in the Hamiltonian. However, in reality, physical systems are also 
subjected to time-dependent noise (this is, for instance, the case of applications to nuclear 
physics [0). In the KR, the influence of a time-dependent noise was discussed for the first 
time in |55|] (see also the detailed investigation in |3"6|), finding that if the strength of the 
noise exceeds some critical value, then it destroys coherent effects of quantum localization 
and pure classical diffusion is recovered. 

The addition of noise to the BRM provides an alternative method to determine the 
diffusion coefficient D from the direct computation of the linear growth of M(b,t). It is 
first interesting to notice that the computation of D cannot be easily performed in practice 
without the presence of noise. The reason is that in the absence of a time-dependent noise, 
corrections to the linear behaviour of M(b, t) arise already at short times, thus preventing 
an accurate determination of the coefficient of the linear growth. In other words, coherent 
effects of quantum localization come into play even on the time scale of classical diffusion 
and this results in a rather smooth transition from classical diffusion to complete relaxation 
(see Fig. 3). 
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In order to estimate the critical value of the noise-strength g 2 r which destroys completely 
quantum coherence and leads to pure classical diffusion, one needs to compare the shift 
of levels induced by the additional noise, with the mean level spacing between operative 
eigenstates. Since the latter turns out to be proportional to 1/b 2 , one can see that if the 
shift AE « g 2 tn is larger than 1/b 2 , localization will be completely destroyed. Accordingly, 
effects of quantum coherence should be observable only when the condition g 2 > 1/b is 
satisfied. Our numerical simulations confirm this estimate: the data in Fig. 4, which refer 
to b = 8, show that D/b 2 rs 1.8 in good agreement with the value found from the fit of the 
scaling dependence given by (|16|). We would like to stress that diffusion due to the noise 
occurs also for g < g cr however, the rate of such diffusion is different from that given by the 
classical diffusion determined in the limit b — > oo (for details see, e.g., f36fl ). 



In any statistical process, the analysis of average quantities provides a limited description 
of the underlying properties. At least the variance should be considered in order to achieve a 
more complete characterization of the phenomenon of interest. In the present case, we shall 
consider the fluctuations of M(b,t) in both the diffusive (t& < £ < £d) and the relaxation 
(t ^> t D ) regime. The relevant quantity to be determined is the size of sample-to-sample 
fluctuations, 



where the brackets (...) denote the average over different realizations of the Hamiltonian 
H nm . A meaningful way to present the numerical data is through the relative amplitude 
H = AM/M of the fluctuations. Its scaling properties with b are discussed in the two 
following subsections. 



IV. FLUCTUATIONS OF PACKET WIDTH 




(19) 
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A. Diffusive time scale 



As the diffusive time scale is relatively short, we have been able to determine \i for quite 
large b- values (namely, b < 45), so that finite band-size corrections should be definitely 
negligible. In our numerical experiments, we have integrated Eq. (|6|) up to a time t m = b 2 /10 
for more than 1000 realizations of the disorder. The results are reported in Fig. 5 under the 
scaling assumption 

/i(M) -b^Jiit/b 2 ) . (20) 

By performing a least square fit of log/i versus log b at fixed time t = b 2 /10, well inside the 
diffusive regime, we have estimated 77 which turns out to be approximately 0.7 (see the inset 
in Fig. 5). One should notice the substantial disagreement with the value of the exponent 
recently obtained in the KR |37 ], rj « 1.0. This scaling parameter has been conjectured 



to be related to the mesoscopic fluctuations of the diffusion coefficient p7fl - However, the 
connection has not been entirely clarified. 

B. Relaxation time scale 

The next important issue concerns fluctuations around the so-called steady state distri- 
bution of the wave packet in the asymptotic regime t^> tp. If one assumes that the steady 
state distribution of c n (t — > 00) is characterized by an ergodic spread of the packet over 
some finite size N s , and if the components c n are statistically independent, then 

AM 1 

a = — — « -= . 21 

M \[W S 

Since the components c n are directly related to the amplitudes of eigenstates and the latter 
are expected to be random on the scale of their localization length, one can conclude that 
N s ~ (^oo(fr)) ~ b 2 , i.e. fi w 1/6. Surprisingly, our numerical data in the range b = 4 4- 12 
indicate that the scaling has a definitely different form, namely 

^(t/b 2 ) « A (22) 



with 5 ~ 0.7 (see Fig. 6). A detailed study shows that when the value of 5 is varied by 
±0.1, the superposition of the various curves on the plateau gets appreciably worse. As 
an additional check, we have performed a least square fitting of fi versus b after averaging 
the curves over times t > 20b 2 . The fit presented in the inset of Fig. 6, confirms the value 
5 ~ 0.7 ± 0.01 . One should note that this value is in perfect agreement with the result 
found on the diffusive time scale for the factor rj, although there is no reason a priori to 
expect such an equality over time scales where different physical mechanisms control the 
packet evolution. We would like to stress that the anomalous scaling described by Eq. (p2|) , 
is in close agreement with the numerical data for the KR, where it has been observed that 
ji ~ b~ s with an anomalous exponent 5 ~ 0.6 [14]. Note that in the KR, larger values of the 
effective parameter b have been reached. 

In order to better disentangle the question of fluctuations, we have investigated also the 
temporal behavior of M(b, t). The main motivation for this study is the comparison between 
sample-to-sample and temporal fluctuations for a typical realization of the disorder. In 
practice, we have integrated the Schrodinger equation for a time up to t = 1256 2 , discarding 
an initial transient time t < to = 20b 2 (which is sufficiently long for M(t) to saturate). The 
Fourier power spectrum |[/(co>)| 2 of u(t) (let us recall that M(t) = (u(t))) signal has been 
then averaged over more than 150 realizations. The results for b = 4, 6 and 8 are reported 
in Fig. 7 in a doubly logarithmic plot. The best data collapse is obtained by assuming 
that (|[/(c<j)| 2 |) ~ 6 ? with <; = 6.6. At "high" frequencies, (\U(lj)\ 2 ) exhibits a Lorentzian- 
type behaviour, which turns, at low frequencies, into a weak divergence that reveals the 
presence of non-trivial long-time correlations. In fact, by invoking the Wiener-Kintchin 
theorem, the low-frequency tail in the spectrum of (|£7(u;)| 2 ) can be connected with the 
relaxation properties of M(t) towards its asymptotic value M^. More precisely, the power- 
law convergence of the type assumed in Eq. ( |16D implies a power-law divergence as 
cd~ l+ P which is compatible with our data. However, the low-frequency cut-off due to the 
finite time of our simulations prevents drawing a definite statement about the presence 
of a truly power-law divergence. Nevertheless, we can, at least, determine the cross-over 
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frequency u c separating the two temporal regimes, which turns out to be u c ~ 0.02/6 2 . 
Such a frequency is approximately 100 times smaller than the mean spacing between the 
energy levels of the eigenstates which effectively participate to the evolution of the wave 
packet. This observation can be taken as an indirect confirmation of the role played by 
Mott states in the long-time evolution. As it has been already recalled, Mott states have 
quite a specific structure: they appear in pairs and are characterized by two humps a distance 
L apart. This leads to a quasi- degeneracy of the order of AE pa exp(— L/^), where L is 
typically much larger than the localization length l^. Accordingly, over long time scales, a 
few Mott states may dominate the packet dynamics. 

For what concerns the scaling behavior of the spectrum (|£7(u;)| 2 ) with respect to b, the 
estimated exponent q pa 6.6 is in perfect agreement with (p2[). Indeed, the relation between 
the spectral density (or power spectrum) and the variance of the signal M(t), implies 

(AM(t)) 2 = Y,(\U(lu)\ 2 ) = Y,b 6S (\U(iub 2 )\ 2 ) = b 6S (AM(t/b 2 )) 2 . (23) 

By comparing Eq. (|23|) with Eq. fl22|), one obtains again 5 ~ 0.7. 



In Ref. [14] it was conjectured that the above anomalous scaling can be considered as 
an indication of the fractal structure of the quantum steady-state distribution c n (t — > oo). 
More precisely, they argued that the asymptotic shape can be described by an ensemble 
of only N s ~ statistically independent degrees of freedom (N s being the number of 
"channels" where the amplitude of the wavepacket is essentially different from zero). The 
same conjecture can be raised in the present case as well, although a direct check is an 
extremely hard task. 

V. STEADY-STATE DISTRIBUTION 

A. General discussion 

As mentioned above, the localization of all eigenf unctions (see Eq. (|J)) implies that for 
t ^> t D the quantum steady-state f(n,t) = \c n (t)\ 2 fluctuates around an average profile 
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f s (n) = (f(n,t)). As the effective number of eigenstates composing a single wavepacket is 
finite, the average profile does depend on the disorder realization. However, in the limit 
b — ► oo, the number of statistically independent components diverges and sample-to-sample 
fluctuations are expected to vanish. In that limit, temporal and ensemble averages should 
coincide as long as the motion is ergodic. 

On the basis of diagrammatic techniques, many results have been obtained for the steady 
state distribution f s (x) in continuous ID models with white noise potential (here, x denotes 
the position of the electron). In particular, an exponential decay f s (x) ~ exp(— \x\ /4l m ) 
has been predicted for the tails of f s (x) (l m , being the mean free path) |24|]. A subsequent 
more accurate analysis |38| revealed the presence of the prefactor |xp 3 ^ 2 . Both findings are 



included in the global expression derived in |3 

j-w - i£r^ v+ *> • < 24 » 

In fact, the above expression implies that, close to the origin, the spatial dependence is 
purely exponential, 

f s (x) ~ exp(-\x\/l m ) ; x < l m (25) 

while the asymptotic decay is described by 

f a (x) ~ |xr 3/2 exp(-|x|/4/ m ) ; x > 4Z m . (26) 

Therefore, the above two equations reveal that the decay rate S(x) = (In f s (x))' (the prime 
denotes derivative with respect to the argument), changes by a factor 4. One consequence 
of the non purely exponential behavior is that the average size of the saturated packet is 
two times smaller than the asymptotic localization length (\x\) = 2l m . It is interesting 
to notice also that the asymptotic dependence (1261) is similar to that for the conductance 



of ID samples of finite size [|0| in the strong localized regime; in this case x is the ratio 
of the sample size with the localization length. As no analytical results are available for 
quasi- ID systems, in the next section we shall compare our numerical results with the above 
expressions, by fitting the only free parameter l m . 
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In any case, some information on the steady state distribution f s (n) can be obtained 
from the structure of the eigenstates by exploiting the following equality, 

fs(n) = J2\<p m (n )(p m (n)\ 2 , (27) 

m 

where tp m (n) is the n-th component of the eigenstate with energy E m and n is the position 
of the initial 5— like packet. Therefore, determining the asymptotic shape of a wave packet is 
tantamount to determining the average correlation properties of single eigenstates. Although 
no rigorous results are known in this direction, a phenomenological approach allowed to shed 
some light on the closely related KR problem P, p0| , p2| . Because of the analogies between 
BRM and the KR, it is instructive to compare the results for the steady distribution as well. 
A rough estimation of the tails of f s (n) can be obtained in the following way (see |T3|J41 



for details). If one assumes the simple exponential form (f m (n) ~ exp(— |m — n|/ioo)> for the 
m-th eigenstate, Eq. ( pT|) leads to the expression 

/.(„)- .^-fc^l) , (28) 

which, in turn, implies l s = l^, where l s is the localization length of the asymptotic packet 
(defined from the probability amplitude, i.e. taking the square root of Eq. (|2"B|)). However, 
this result is inconsistent with the numerical data for the KR which instead indicate that 



l s ~ 2/oo |yj. To explain the latter result, it was suggested to take into account the large 



fluctuations of the eigenstates around their shape, 



1 / 1 771 — Tl I \ 

fm(n) = exp Vimn) , (29) 



where £ mn is a Gaussian noise with the zero mean and a variance [II], 



< (A£ m „) 2 > = D 3 \m-n\, D s ~ 1//^ . (30) 

By inserting the Ansatz ( p9|) in Eq. (|27|) and averaging over the noise term, it was found 
that 

( \n — m\\ 

(\(Pm{n)\) ~expl I , (31) 
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implying that the linear average yields a different localization length compared to (|28|). 

By repeating the same calculations for the expression (p7[), one obtains that f s ~ exp(\n— 
no\/loc,), which is found to be in agreement with the numerical results for the KR model as 
it implies l s = 21^. 

However, relation (f27|) implies that the average < |0 m (n)| 2 > should be used rather than 
{\4>m(n)\). In such a case, the result is 




which implies that f s ~ exp(\n — nol/2/oo). The correctness of this expression is confirmed 
by the relation l s = 41^ that it implies. 

B. Numerical data 

In order to determine the asymptotic shape of the wave packet, we have followed the 
evolution of an initially 5-like packet for times t > 120ti> The distribution f s (n) has been 
then obtained by averaging over more than 150 realizations for several 6-values in the range 
b = 4 12. The results are reported in Fig. 8 with the by-now-standard scaling hypothesis, 

f s (x) = b 2 f s (n), x = n/b 2 , (33) 

that is once more confirmed by the good data collapse. 

A peculiarity of all our simulations is that f s (no) is larger than the neighboring values 
by approximately a factor 3. The reason of this apparent anomaly can be traced back to the 
specific 5-like shape of the initial packet (that implies Eq. (|27|)) and to the spatial random 
structure of the eigenvectors. Indeed, the latter assumption, together with the observation 
that only a finite number L of channels effectively contribute to the sum in Eq. (p7|), leads 
to 

/.(no) « L ( V \n )) ; f s (n + n ) « L 2 (<p 2 (n)) (<p 2 (n)) . (34) 
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Since it is known that (ip A ) = 3/L, while (if 2 ) = 1/L for the eigenfunctions of matrices 
belonging to the Gaussian Orthogonal Ensemble ||, we obtain that f s (n )/ 'f s (n) = 3 for n 
close to but different from uq. This value of the ratio is in pretty good agreement with the 
above mentioned numerical estimate. 

Moreover, the numerical results reported in Fig. 8 strongly suggest that the decay of 
f s {n) in the vicinity of no is definitely faster than in the tails. Therefore, it is very tempting 
to compare this data with the theoretical dependence derived for ID disordered models 
(Eq. (|24f) ). The best fit of the only free parameter gives / m ~ 0.29; the corresponding curve 
is shown in Fig. 8 (see the solid line). The very good agreement between the numerical 
results and the analytical curve over a broad range of x-values suggests that a properly 
modified theory to include the determination of the mean free path from first principles 
should be able to account for the asymptotic properties of packets in quasi-lD systems as 
well. 

One should also notice that the dependence of the slope on the distance from the center 
n of the packet is an entirely new feature with respect to the analogous problem in the KR 
4"T| , |T3"|] , where no evidence of two distinct regions of localization has been found. The reason 



of this discrepancy is not clear: on the one hand it is possible that this reflects an actual 
difference between the two models, on the other hand, it is possible that the accuracy of 



the numerical data for the KR |4T| is not enough to reveal this peculiarity in the associated 
profile f s (n). 

A further difference is the variation of the localization length of the eigenstates with 
the energy in BRM (see Eq. (|5D). Far from the center of the packet, we expect that the 
decay is dominated by the longest localization length Zoo(0) = 2/3 (in b 2 units). On the 
other hand, from Eq. (|26|), we find that l s = 8l m ~ 2.32, which is only slightly smaller 
than 4^(0) ~ 2.66. Accordingly, the equality l s = 4/^ found in the phenomenological 
theory for KR and explained by invoking the presence of strong fluctuations of the individual 



eigenstates ET| appears to hold also in the present case. The small deviation is presumably to 



be attributed to the not yet vanishing contributions of more localized eigenstates. Instead, 
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if we average over all energies, we obtain (l^E)) = 0.5 (in b 2 units) which results to a 
localization length l s = 4(Z 0O (£ 1 )) = 2 for the total wavefunction. Moreover, it is interesting 
to notice that a direct determination of l s by fitting the profiles reported in Fig. 8 with a 
pure exponential law, yields l s ~ 2: this means that the multiplicative correction l/|x| 3 / 2 
in expression (^Bj) is essential for a correct estimate of l s , if the range of n/b 2 is not large 
enough. 

A further more direct confirmation of the presence of strong fluctuations in the various 
eigenstates is obtained by determining the localization length lj® of the asymptotic packet 
from the logarithmic average ((ln(/ s ))) of the single packets. The resulting profile, reported 
in Fig. 9, yields ~ 1.3 to be compared with the value l s ~ 2.32, obtained from the 
arithmetic average. Interestingly enough, the ratio between the two lengths and l s for 
the shape of saturated wavepacket, is approximately equal to 2 as known for the single 
eigenstates (see the previous subsection). 

C. Fluctuations of the steady-state 

In this section we investigate directly the nature of wave-packet fluctuations in the asymp- 
totic regime. This enables us to test the correctness of a conjecture relative to the single 
eigenvectors raised in the context of the KR. Indeed, fluctuations are very important in that 
they allow explaining the difference between the decay rate of the wave packet and that of 
the eigenvectors. 

Having in mind Eq. (p9j), we computed the logarithm of f s (n) and studied its variance 
at a distance An = 6b 2 from the center of the packet (averaging also over a small window of 
5 neighboring sites, under the assumption that the fluctuations are nearly constant in such 
interval). As a result, we have found that the distribution function of y = ln(/ s ) is, with a 
good accuracy, a Gaussian (see for instance the histogram reported in Fig. 10 which refers 
to the case 6 = 5 and is the result of 5,000 simulations with independent realizations of the 
disorder). This represents a first confirmation that the hypotheses made in the KR can be 
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profitably carried over to the present model. 

A more complete information is obtained by studying the fluctuations of y for different 
values of An. More precisely, we have computed the variance 

a 2 (n) = ((ln/ s (n)) 2 )-(ln/,(n)) 2 (35) 

in the steady-state regime. The results for the cases 6 = 5^-8 are reported in Fig. 11 with 
the scaling assumption a 2 (n) = b 2 a 2 (x), where x = An/6 2 . The data shows that for large 
x, the variance grows linearly with x, indicating that the logarithm of the profile diffuses 
around the average value. This is a further confirmation of the validity of a relation of the 
type (j29D for the wave packet and it is strengthened by the observation that the estimated 
slope is approximately equal to 1 in scaled units. 

This behaviour is also analogous to what happens in ID disordered models, where it is 
assumed that the logarithm of the absolute value of the Green function In \G(m,n; E)\ has 
a Gaussian probability distribution with the mean equal to — \m — nl/l^ and the variance 
equal to \m — nl/l^ . 

Another interesting observation concerns the central part of the packet (i.e., |x| < 1), 
where the variance a 2 is almost constant, indicating that the amplitude values are essentially 
independent of one another (see also the inset in Fig. 11). Accordingly, all numerical findings 
do confirm the conjectures that have been so far utilized to present all the features of 
wavepacket diffusion in a coherent manner. Unfortunately, so far there are no analytical 
results concerning the structure of the eigenstates in the middle of their localization region, 
even in the well-studied ID Anderson case. 



VI. CONCLUSIONS 

In the present paper we have studied the evolution properties of wavepackets in quasi 
ID disordered media described by tight-binding Hamiltonians with long-range random in- 
teractions. We have found that the wavepacket: (i) first spreads ballistically, over a time 
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scale of the order t ~ 1/5 ' 5 , which becomes negligible in the limit b — > oo; (ii) exhibits a 
diffusive behaviour, for times of the order t ~ b 2 \ (iii) finally, for times larger than to > b 2 , 
stops spreading remaining asymptotically localized. 

The scaling properties of the spread of the packet are different in the ballistic (see 
Eq. ([TOD) and diffusive (see Eq. flT2"|)) regime. Beyond the ballistic regime, we propose 
the heuristic formula (|16|) to effectively describe both the diffusive spreading and the even- 
tual saturation. The most interesting feature of Eq. (|i~6|) is the prediction of a power-law 
convergence of the wavepacket width M to its asymptotic value, the deviation going to zero 
as 1/t 13 with (3 ~ 0.8. The parameter (3 accounts for the repulsion between those eigenstates 
which are effectively excited by the initial wavepacket. More precisely, the slow convergence 
of M is attributed to effect of Mott states, that are expected to give eventually rise to a 
logarithmic time dependence, that cannot be numerically observed. 

The "short" time diffusive regime has been investigated by computing dM/dt when a 
small amount of time- dependent disorder was superimposed to the quenched disorder. As 
a result we found that the diffusion coefficient in dimensionless units is approximately four 
times larger than the localization length, after averaging over the energy dependence. The 
presence of noise with strength g 2 of the order of a critical value g 2 r ~ 1/6, destroys quantum 
coherence and recovers classical diffusion. 

Another issue addressed in this paper concerns the fluctuations of the size of the packet 
over different time scales. We found that for both diffusive and saturation time scales, 
the relative amplitude of the fluctuations scales as fi = AM/M ~ l/b v where 77 ~ 0.7. 
We confirmed this law by taking the Fourier power spectrum |[/(c<j)| 2 of sample to sample 
fluctuations of M. This result, besides confirming the anomalous scaling behaviour provides 
evidence for the ergodicity of the evolution. 

The study of the Fourier power spectrum shows that the tails of |£7(u;)| 2 have a 
Lorentzian-type form. At small frequencies, a weak singular behaviour for |£7(u;)| 2 has been 
detected; this result is in agreement with the power-law convergence of M to its asymptotic 
value. Furthermore, we have been able to extract the cross-over frequency u c , where devia- 
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tions from 1/uo 2 start in |£/(u;)| 2 . Such a frequency is two orders of magnitude smaller that 
the mean average separation between the energy levels which participate to the evolution 
of the wavepacket. This is a further confirmation of the participation of Mott states to the 
long-time evolution. 

For what concerns the asymptotic shape of the wavepackets, we found that the scaling 
law (^) is well verified already for relatively large b values. Moreover, the analytic expression 
derived in the context of the ID Anderson model (see Eq. (|24f) ) reproduces pretty well the 
shape of the average profile, upon fitting a single parameter. A further interesting result 
of our numerical analysis concerns the difference by a factor 4 between the decay rate of 
the asymptotic profile around the origin and that along the tails. Moreover, the difference 
between the localization length estimated from the average of the logarithm of the profile and 
that obtained from the linear average of the profile confirms the relevance of fluctuations. 
This observation is reinforced by the asymptotic linear growth of the variance a 2 of the 
logarithm of the profile versus the distance from the center of the packet. 
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Figure Captions 



Fig. 1 Scaling of M(b, t) vs. time for the ballistic time scale. The reported values of b 
are b = 20, 25, 30, 35, 40, 45. The inset shows the lost of scaling for t > i&. 

Fig. 2 An example of wave-packet diffusion beyond the ballistic time scale for b = 12. 
The same curve is reported (in a doubly logarithmic plot) in the inset, magnifying the early 
stages of diffusion, to testify the linear behaviour. 

Fig. 3 Mean square displacement M for b = 4, 5, 6, 7, 8, 10, 12. The smooth curve 
corresponds to the phenomenological expression (16). In the inset the same quantities are 
shown for shorter times. 

Fig. 4 Diffusion of packets with additional noise for 6 = 8. The diffusive constant is 
determined from simulations performed for different noise strengths g 2 . 

Fig. 5 Relative fluctuations on the diffusive time scale. The data is scaled with respect 
to the law (20). In the inset, a least square fit for the same quantity and different b- values 
is shown; the time is fixed, t = 0.16 2 . 

Fig. 6 Relative fluctuations in the saturation regime. The values of b are 

4, 5, 6, 7, 8, 10, 12. In the inset, a least square fitting is shown for the same quantity which 
now is averaged over time from t = 20b 2 up to t = 1206 2 . 

Fig. 7 The Fourier spectrum (\U(u))\ 2 ) for the cases b = 4, 6, 8 and for times t > 20b 2 . 

Fig. 8 Asymptotic average profile of the wavepacket for b = 4 4- 12 after rescaling. The 
inset shows the behavior near to the maximum; smooth lines follow from the theoretical 
expression (|24|) . 

Fig. 9 Logarithm of the steady-state f n for b = 5. 

Fig. 10 Distribution of log f n for 6 = 5. 

Fig. 11 Variance of log/ n in the steady-state after rescaling the for b = 

5, 6, 7, 8, 10. The straight line is the fit for large values of x. The inset shows the behavior 
close to the maximum. 
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